Send nudis: An assessment of nudibranch diversity in Sodwana Bay, South Africa

Abstract Climate change is posing unprecedented pressure onto marine ecosystems worldwide. This makes it imperative to monitor the effects that are being experienced in these environments. Nudibranchs are benthic marine organisms that possess characteristics that have the potential to act as indicators of change within ecosystems such as coral reefs. Therefore, these species have the ability to provide valuable information on fine‐scale changes in environmental conditions. It is thus essential for studies, such as this, to establish baseline analyses from which changes within nudibranch populations can be examined in order to investigate their ability to act as bioindicators. Recommendations can also be made for future sampling procedures through investigating environmental and experimental parameters that influence nudibranch communities. Nudibranch populations were sampled on Two‐Mile Reef in Sodwana Bay, South Africa, through SCUBA where individuals were photographed and later identified. Data were collected within a sample‐based dataset, as well as by citizen scientists within an incidence‐based dataset. Across both datasets, a total of 85 species were identified. Nudibranch populations showed high levels of diversity within an uneven, unstable community. Citizen scientist data provided imperative information to the baseline assessment and, therefore, the inclusion of these data increased the robustness of this study. Environmental and experimental variables investigated did not influence the outcomes of this study and should therefore not be heavily focused on in designing future experiments. Future monitoring studies should continue to record oceanic pH in order to detect any possible changes due to ocean acidification. It is recommended that sampling events should be increased in order to capture all species present in these localities. These events should also encompass an extended temporal scale in order to cover a larger temperature range. Research on bioindicators is essential within today's rapidly changing climate, mainly due to human activities, particularly within an extremely vulnerable habitats such as coral reefs.


| INTRODUC TI ON
Resource availability of all ecological communities is chronically altered by drivers of global change (Koerner et al., 2015).
Resources such as nitrogen deposition, deviations in precipitation, and increasing carbon dioxide levels are being impacted (Koerner et al., 2015).Anthropogenic activities are accelerating the rates of global change at an alarming speed by intensifying the effect of these drivers (Ingole, 2021).It is anticipated that these alterations in resources will have large impacts on ecosystem functioning through transforming community structure and composition (Koerner et al., 2015).Marine ecosystems form intricate food webs and are thus characterised by complex biological interactions that directly or indirectly impact the success of species through the performance of another (Doney et al., 2012).The main physiological effects that rising atmospheric carbon dioxide levels and global climate change will have on marine environments will include concurrent increases in temperature, changes in circulation, increased stratification, changes in nutrient input, ocean acidification, and decreased oxygen content (Brierley & Kingsford, 2009;Doney et al., 2012).These changes are expected to affect the natural composition of marine populations by causing phenological disruptions, shifts in abundance, as well as changing the spatial organisation of organisms such as distribution and dispersal (Brierley & Kingsford, 2009;Doney et al., 2012).An effective method for assessing the effects that ecosystems experience due to climate change can be achieved through the use of bioindicator species (Cooper et al., 2009).
Bioindicators are organisms that reflect the health or quality of an ecosystem (Jackson et al., 2000;McGeoch, 1998;Siddig et al., 2016).They are therefore frequently used by natural resource managers for assessing the well-being of ecosystems, as well as providing early warning signs of environmental deterioration (Jackson et al., 2000;McGeoch, 1998;Siddig et al., 2016).
The presence of bioindicator species can convey important information about biological, chemical, or physical changes occurring in the environment and can therefore indicate major stressors experienced by ecosystems (Niemi et al., 2004;Niemi & McDonald, 2004).The use of bioindicators is becoming increasingly prevalent with natural resource management agencies as they provide an easy, rapid, and inexpensive method in depicting environmental degradation on a larger scale (Cooper et al., 2009;Navarro-Barranco et al., 2020).It is anticipated that marine bioindicator organisms will be crucial in the conservation and management of coral reefs worldwide by highlighting changes within the physical environment (Dale & Beyeler, 2001;Navarro-Barranco et al., 2020;Niemi et al., 2004).Due to their unique life history traits and ability to indicate change, this study sets out to suggest the use of nudibranchs as essential bioindicator organisms to understand the effects of climate change on coral reef ecosystems (Goddard et al., 2011;Nimbs et al., 2015).
Nudibranchs are a diverse group of brightly coloured marine gastropod molluscs that shed their shells in the larval stage of development and thereafter remain shell-less (Dean & Princep, 2017;Dumdei et al., 1989).There are over 4700 known species of nudibranchs (Gofas, 2021) which inhabit a diverse array of habitats.
As they are essentially blind, nudibranchs perceive their environment through two specialised rhinophores that are chemosensory structures (Arey, 1918;Murphy & Hadfield, 1997).The rhinophores are used to interpret chemical cues transported by water currents in order to locate food or potential mates (Arey, 1918;Wertz et al., 2006).Due to their reliance on chemoreception, nudibranchs are sensitive to changes in ocean water parameters that are affected by global warming, such as changing pH associated with ocean acidification (Albright et al., 2018;Kurnianda et al., 2020;Seroy & Grünbaum, 2018).Nudibranchs are also susceptible to environmental changes that affect the substrate on which they depend on for feeding and reproduction (Biermann et al., 1992;Burn, 2015;Johnson & Gosliner, 2012).
Species in this order are carnivorous benthic organisms that prey on coral, sponges, hydroids, sea anemones, fish eggs, barnacles, as well as other nudibranchs that they can locate on the sea floor (Kara et al., 2018).With only a few known exceptions, nudibranchs are generally seen to be hermaphrodites (Burn, 2015;Sekizawa et al., 2013), with both male and female apertures located within a common chamber called an atrium on the anterior right side of their bodies (Burn, 2015).During mating, individuals align so that there is contact between their reproductive apertures (Burn, 2015).
Nudibranchs then deposit spirally coiled gelatinous egg ribbons onto hard surfaces (Biermann et al., 1992;Burn, 2015).The adult nudibranchs die once the egg-laying process is complete; thus, the life-cycle is short-lived (Burn, 2015;Costello, 1938).Smaller species have a life span of 4-6 weeks, whereas some larger species may live approximately 2-3 years (Burn, 2015).Their short life span allows them to respond rapidly to changing environmental conditions (Kurnianda et al., 2020).Overall, due to their life history traits, reliance on chemoreception, and their ability to rapidly adapt, nudibranchs make for key bioindicator species of reef health, as well as environmental changes (Goddard et al., 2011;Nimbs et al., 2015).These traits make long-term data on diversity, abundance, and biogeography of nudibranchs essential for understanding the effects of climate change on marine ecosystems.

Biodiversity ecology
Many studies to date have been conducted almost entirely on the bioactive metabolites and taxonomy of nudibranchs, whereas only a few studies have focused on diversity and abundance (Sabdono et al., 2021).None of these limited studies have been conducted in Sodwana Bay, and therefore, this study is the first of its kind.The importance of this study is that it lays the foundation for future studies to investigate nudibranchs as bioindicators of climate change, as well as the effect that climate change will have on the coral reefs of Sodwana Bay.
In order to establish nudibranchs as bioindicators in future research, the main aim of this study was to provide a baseline assessment of nudibranch diversity on Two-Mile Reef, a marine protected area, located in Sodwana Bay, South Africa.This was assessed through an underwater visual census of nudibranchs at three dive sites.Citizen scientist data of nudibranch observations were analysed in parallel with scientific findings in order to increase the robustness of the census.Additionally, the influence of environmental and experimental parameters on nudibranch observations was investigated.This was done to determine if sampling procedures are affected by these factors.Finally, the results of this study will be used to recommend the most efficient sampling methods to assess nudibranch diversity for long-term monitoring programs.Scientific findings will be logged into a sample-based (SB) dataset and citizen scientist data into a separate incidence-based (IB) dataset.Within the SB dataset, it is predicted that diversity among the three study sites will be similar, as indicated by the Shannon-Weaver index.It is also anticipated that the uniformity index will show that an unstable community is present at all three sites, due to differences in abundance among species.Chromodoris hamiltoni Rudman, 1977, and Chromodoris celinae Tibiriçá, 2019, are foreseen to be the most common species at all the sites and will therefore affect the Simpson Dominance index.Other species expected to be observed include

| Location and timing
Sample-based (SB) observation data were collected in the marine protected area of Sodwana Bay, South Africa across three different dive sites.All three of the study sites are located on Two-Mile Reef.These dive sites were Chains, Anton's, and Bikini reef with the maximum depth for each study site at 17.00, 15.00, and 21.00 m, respectively (Figure 1).Data collection was conducted

| Sample-based data collection
Underwater visual census were utilised to obtain data on nudibranchs.This involved a team of SCUBA divers observing nudibranchs along a designated route set out across each dive site (Figure 2).Each study site was sampled five times over the period of 2 weeks.Each nudibranch observed was photographed using a Nikon Coolpix w300 underwater camera.Species were then identified using Nudibranch & Sea Slug Identification -Indo-Pacific 2nd Edition (Gosliner et al., 2015), A Guide to the Sea Slugs of the Maputaland Coast (Strö ̈mvoll & Jones, 2019), The Reef Guide (King, 2014), The Sea Slug Forum online database (Australian Museum, 2010), and the Nudibranch ID Indo Pacific application (Cobb, 2022).Journal articles were also utilised for identifying species within the Chromodoris (Tibiriçá et al., 2019) and Flabellina (Ekimova et al., 2022) genera.Species were logged into a dataset for each study site.Various parameter measurements were collected in situ for each of the dive sites using a Thermo Fisher Scientific handheld conductivity metre and a Suunto Zoop dive computer.Environmental parameters collected included temperature (°C) and acidity (pH).A purposive sampling method was used to obtain water parameters, such as visibility, and the strength of the surge and current.Experimental parameters such as dive time, maximum depth, visibility, and number of observers were also recorded for each sampling event.

| Incidence-based data collection
Observations from citizen scientists in Sodwana Bay were obtained through the Deep Blue Conservation organisation.Nudibranchs were also photographed with a Nikon Coolpix w300 underwater camera.Identification of species were carried out using the same resources as utilised in the SB dataset.Environmental parameters collected by citizen scientists included temperature (°C), strength of the surge, and the strength of the current experienced during each dive.Experimental parameters collected included dive time, maximum depth, visibility, and number of observers.

| R analysis
Due to the different sampling methods employed for the SB data and IB data, they were treated as separate datasets in all the analyses.R version 4.1.2(R Core Team, 2021) was used for all statistical analyses in this study.Generalised linear models were used to test for significant differences between species richness and abundance between sample sites in the SB dataset.A Venn diagram was produced to show common or unique species belonging to the two datasets.Species accumulation curves were generated for each dataset.Sample-based data were used to create a species accumulation curve based on sampling events.The number of individuals was used to produce the species accumulation curve based on IB data.A rank abundance curve of SB data was produced for all species observed across the three dive sites.A second rank abundance curve of IB data was also produced.The Shannon-Weaver index as well as the Simpson dominance index were calculated in RStudio using the "diversity ()" function.These values were obtained for each sampling event on each of the dive sites and an average was then obtained for the SB dataset only.The uniformity index values were calculated manually with the formula below using the average Shannon-Weaver index for each dive site.Generalised linear models were constructed using the environmental and experimental parameters as predictor variables, with nudibranch abundance or species richness as the response variable.An ANOVA was then performed to test for a significant effect of the predictor variables on the response.Significance of both environmental and experimental parameters on species composition was also obtained using a permanova test.

| Shannon-Weaver index
The Shannon-Weaver index (H′) was used as an index of diversity in this study.The H′ was calculated for each sample site in order to provide a baseline diversity level for Chains, Anton's, and Bikini in the SB dataset (Clarke & Warwick, 2001).The formula used was (Shannon & Weaver, 1998): where Pi = ni/N, with ni as the number of individuals of each species (ith species) with N as the total number of individuals in the sample (Paulangan et al., 2021).The results obtained for the Shannon-Weaver index was compared to the criteria for assessing the results (Table 1).
F I G U R E 2 A map of each dive site for the sample-based dataset; Chains, Anton's, and Bikini Reef.The red broken line indicates the route that was followed each time the dive site was sampled.The start and end point are depicted on the route as well as an arrow to indicate the direction of travel.

TA B L E 1
The criteria for interpreting the Shannon-Weaver index (H′) of diversity utilised in this study (Paulangan et al., 2021;Shannon & Weaver, 1998).

| Uniformity index
The uniformity index (E) describes the number of individuals between species within a community (Ulfah et al., 2019).If the individuals between species are evenly distributed then the more balanced the ecosystem will be (Ulfah et al., 2019).The formula used is (Odum, 1967): where H′ is the Shannon-Weaver index and H max = ln(S), where S is the total number of species found (Krebs, 1999).The uniformity index falls within a range 0-1 and follows the proceeding criteria (Table 2).
A low uniformity index shows an uneven distribution of species within the community where one species tends to dominate (Ulfah et al., 2019).A high uniformity index shows evenness between the number of individuals of each species that inhabit a community (Ulfah et al., 2019).

| Simpson dominance index
The Simpson dominance index (C) was used to assess the evenness of species composition within a community (Paulangan et al., 2021).
The formula used (Odum, 1967) is as follows: where Pi = ni/N, with ni as the number of individuals of each species (ith species) and N as the total number of individuals in the sample (Odum, 1967).Much like the uniformity index, the dominance values also range between 0-1 and are assessed according to the following criteria (Table 3).unique to the SB dataset and 22 to the IB dataset, bringing the total identified species to 85 (Figure 3).

| RE SULTS
The species accumulation curve generated from the SB dataset (Figure 4a) showed incomplete sampling.An asymptote has not been reached within this dataset showing that not all species that are present were sampled.This curve indicates that the total species observed would increase until reaching an asymptote if additional sampling events were undertaken.Similarly, the species accumulation curve of the IB dataset also depicts incomplete sampling (Figure 4b).
However, the curve is closer to reaching an asymptote when compared to that of the SB dataset (Figure 4a).
The rank abundance curve produced from the SB dataset (Figure 5a) indicates that there are levels of unevenness in the community.Higher ranked species occur at greater abundances and therefore dominate.There is a large disparity between the abundance of the highest ranked species and those that proceed them.
The lowest ranked species occur at similar abundances as the curve remains approximately horizontal as rank decreases.A very similar rank abundance curve was generated from the IB dataset (Figure 5b).
This curve also depicts unevenness within the community and a large disproportion in abundance between higher and lower ranked species.
The diversity indices were calculated for the SB dataset only.
The Shannon-Weaver index ranged between 1.56 and 2.80 across all three dive sites.The average Shannon-Weaver index for Chains, Anton's, and Bikini was calculated to be 2.13, 2.12, and 2.23 respectively.The overall average for the Shannon-Weaver index for all three dive sites was 2.16.The Simpson dominance index ranged between 0.77 and 0.93 across all three dive sites.
The average Simpson dominance index for each of the dive sites Pi 2

TA B L E 2
The criteria for interpreting the uniformity index (E) of community stability utilised in this study (Krebs, 1999).The pH ranged from 7.96-8.17across all three dive sites.The overall average pH was recorded to be 8.07.The surge was never strong and varied between none, low, and medium.The current was only seen to be slight or no current present.Maximum depth was recorded to be within the range of 22.60-15.40m with an average maximum depth of 18.52 m.The longest dive time was recorded to be 75.00min compared to the shortest dive time of 42.00 min.

Stability
The average dive time spent across all dive sites was 58.13 min.The water temperature was not tested for significance on nudibranch assemblages since it was constant at 24°C throughout the study except on the fifth dive on Bikini when it decreased by 1°C.This trend was also seen to be true for visibility which remained constant at approximately 20.00 m.The generalised linear model testing for a significant effect of environmental and experimental parameters on nudibranch abundance, showed that none of the predictor variables significantly affected the response (p > .05)(Table 4).The second generalised linear model testing for a significant effect of environmental and experimental parameters on nudibranch species richness, again displayed that none of the predictor variables significantly affected the response (p > .05)(Table 5).The permanova test showed that no water or dive parameters had a significant effect on nudibranch composition (p > .05)(Table 6).Overall, none of the water or dive parameters assessed in this study significantly affected the nudibranch assemblages sampled.

| DISCUSS ION
This study set out to quantify nudibranch communities within Sodwana Bay, South Africa, to provide a baseline assessment for future studies investigating the use of nudibranchs as bioindicators of climate change.Since this was the first study of its kind, another important aim was to investigate the environmental and experimental

TA B L E 3
The criteria for interpreting the Simpson Dominance index (C) of species evenness within a community utilised in this study (Odum, 1967).there is an uneven distribution of species (Appendix 5).The species that dominated within the SB dataset were H. wasinensis and G. bonwanga.However, in addition to H. wasinensis, C. hamiltoni was highly predominant in the IB dataset.This uneven spread of species is further depicted by the rank abundance curves (Figure 5).These curves show the disparity in abundances of species at higher and lower ranks.Uneven distribution of species, with levels of dominance, is further shown using the uniformity index.This index depicts an unstable community which cooperates with conclusions drawn from the Simpson dominance index and rank abundance curves (Figure 5, Appendix 5).

Simpson dominance index (C) values
Corresponding results from a study conducted in southern Mozambique, 200 km north of Sodwana Bay, found H. wasinensis to be the most common species (Tibiricá et al., 2018).Chromodoris hamiltoni is also known to be a proliferate species found along subtidal reefs located on the east coast of Africa (McPhail & Davies-Coleman, 1997).On the other hand, G. bonwanga is a newly described species (Matsuda & Gosliner, 2018).Thus, there is no data available regarding the abundance or possible reasons for the dominance of this species along the southern African coastline.
Therefore, it is important that aspects like breeding patterns and diet be thoroughly investigated in future studies in order to address  this gap in our understanding of these observations.The diversity indices can be utilised as baseline values for future studies to track nudibranch community level changes.The method of repeated sampling along designated routes was seen to be effective in quantifying nudibranch communities.The success of this was shown through no significant difference between nudibranch abundance or species richness between sample sites.
The average pH recorded was 8.07 which is in line with the current global average of 8.10 (Fassbender et al., 2021).This recorded value can be used to detect minor changes in pH within the Sodwana Bay marine protected area.Oceanic pH levels are expected to decrease in the near future due to the increasing levels of atmospheric carbon dioxide, resulting in an increased level of absorption by the ocean (Feely et al., 2009;Orr et al., 2005).Preliminary readings, as done here, become extremely important in monitoring the effects of ocean acidification on marine ecosystems.
Research has shown the potential of nudibranchs to act as bioindicators of climate change, as well as reef health.This is due to the influence of climate change on nudibranch abundance through the alteration of environmental variables.Total nudibranch abundance is highly correlated with ocean water temperature (Goddard et al., 2011(Goddard et al., , 2013;;Schultz et al., 2011;Wilson et al., 2016) which is shown through population fluctuations in accordance to El Niño and La Niña events (Goddard et al., 2011).During periods of warming, such as those associated with El Niño events and the warm phase of the Pacific Decadal Oscillation, there is a strong positive correlation in nudibranch abundance (Goddard et al., 2011).These climatic events strongly impact sea surface temperature as well as sea surface height, which ultimately effects abundance (Goddard et al., 2011).Contrary to this, abundance is negatively correlated to periods of cooling, such as La Niña events, the cold phase of the Pacific Decadal Oscillation, heightened coastal upwelling, and the positive phase of the North Pacific Gyre Oscillation (Goddard et al., 2011).This dependency on ocean water temperature is precipitating range shifts along with regional declines of certain nudibranch species due to climate change resulting in warming sea temperatures (Goddard et al., 2011;Nimbs et al., 2015).Range shifts Potential draw backs of this study can be attributed to the cryptic nature of nudibranch species (Churchill et al., 2014;Epstein et al., 2019;Fritts-Penniman et al., 2019;Matsuda & Gosliner, 2018).
Due to this, there is potential for cryptic species to be nested within those known to science.Without taxonomic evaluation through genetic assessment or analysis of radular morphology and reproductive systems, the identification of these cryptic species is highly unlikely.
An example of this can be seen in juveniles of certain species of the Chromodoris genus, which are very difficult to distinguish from each other.During early life stages these species do not possess the distinguishing characteristics that the adult individuals show.This was a challenge encountered by citizen scientists that were unable to identify to which species a Chromodoris juvenile should be assigned to.Citizen scientists were also unable to identify species belonging to the Dermatobranchus genus as it is dominated by cryptic species that have similar external anatomy, and can predominantly be distinguished using internal anatomy (Gosliner & Fahey, 2011).With this being said, nudibranch individuals in the SB dataset were identified through a rigorous process using numerous resources such as, identification guides, online databases, applications, expert consultation, and published journal articles, in order to increase the likelihood of correct classification.
Environmental and experimental factors in this study were not seen to influence nudibranch abundance, species richness, or composition (Tables 4-6).These factors included strength of the surge and current, depth, dive time, or number of observers.It can therefore be suggested that these variables should not influence future sampling procedures and do not need to be factored into experimental designs.However, there were methodical limitations present in this study that can be used to enhance future census.Firstly, sampling occasions within the SB dataset were not sufficient to accurately quantify nudibranch populations (Figure 4).Consequently, increasing these events is strongly advised in order to sample all species present.Secondly, temperature remained constant during the sampling period and could not be incorporated into the results.It is suggested that multiple sampling occasions should be completed throughout the year in order to obtain an adequate temperature range.This extension of sampling procedures could also be beneficial as it could encompass more variation within surge and current that was not obtained here.Through the use of photography, distinguishing morphological characteristics could be identified for accurate classification of species that would otherwise be difficult to recall from memory.Photographing individuals is a useful tool and should therefore be continued in future studies.
Data gathered by citizen science has extensively been used by scientists, policy makers, and the general public (Miller-Rushing et al., 2012;Silvertown, 2009).With the aid of citizen scientists, ecological research data can be collected at unprecedented spatial and temporal scales (Dickinson et al., 2012).This is possible through the assembly of historical data, and a large, distributed team of knowledgeable observers (Dickinson et al., 2012).Gathered data can then be employed in critical baseline studies to either respond to crises or identify threats to wildlife and/or people (Dickinson et al., 2012).
Along with this, citizen science encourages the public to become involved in scientific research (Kobori et al., 2016).By integrating the public into such research, education on important matters such as climate change and environmental stewardship is enhanced (Marshall et al., 2012).This study is evident of the important contributions of citizen scientists.A total of 22 new species were added to the baseline assessment through the incorporation of the IB dataset.
Inclusion of this dataset also enabled the investigation of a broader temperature range on nudibranch seasonality.Additionally, volunteers forming part of Deep Blue Conservation were exposed to important concepts such as climate change and using bioindicators to monitor variation in environmental conditions.The incorporation of citizen scientists has been beneficial in this study and should continue to be made use of in future research.
Studies similar to this one are extremely important for laying down a foundation for building, and extending monitoring programs in the future.Within marine environments, these baseline studies are imperative due to degradation being more pervasive than land (Knowlton & Jackson, 2008).Baseline studies are particularly scarce for marine benthic ecosystems such as coral reefs (Knowlton & Jackson, 2008).This lack of information is alarming as these ecosystems are among the most highly threatened by global climate change (Bellwood et al., 2004;Hoegh-Guldberg et al., 2007).To date, many of the ecological studies of coral reefs have been conducted on systems that are already heavily degraded by human impacts (Knowlton & Jackson, 2008).These studies do not provide accurate baseline assessments as they quantify already degraded environments and are therefore not a precise representation of ecosystem functioning (Knowlton & Jackson, 2008).Since Sodwana Bay was declared a protected area in the 1950s, it is a pristine coral reef system (Schleyer & Tomalin, 2000) on which an accurate baseline assessment can be conducted.With nudibranchs at its centre, this study can therefore provide an explicit baseline analysis that can be used to precisely quantify future impacts of climate change on this ecosystem.
Although nudibranchs do not possess all the typical characteristics of bioindicator organisms such as well-studied taxonomy and ecology (Holt & Miller, 2011), their life history traits, reliance on chemoreception, and their ability to rapidly adapt deem them worthwhile to justify future research to investigate their bioindicator potential.

| CON CLUS ION
Nudibranch populations on Two-Mile Reef can be characterised as having a high diversity, uneven distribution of species, and unstable composition.A total of 85 species were identified in both the

ACK N OWLED G M ENTS
We would like to thank Deep Blue Conservation for sharing their data and providing extensive knowledge on the research topic, along with all those involved in collecting these data.We thank Pisces Diving Sodwana Bay for accommodating our research requirements and divemasters for valuable information on the reef systems.We thank Gary Kobb for assisting with species identification.We thank iSimangaliso Wetland Park and Ezemvelo KZN Wildlife for supplying research permits and allowing data collection.We thank Professor Mark Robertson (University of Pretoria) for aiding with data analysis.
Finally, we thank the editors and reviewers for their contributions to previously submitted drafts of this manuscript.

FU N D I N G I N FO R M ATI O N
None.

CO N FLI C T O F I NTE R E S T S TATE M E NT
There is no conflict of interest.

A PPE N D I X 2
Photographs of all species observed within the SB dataset arranged in alphabetical order.Species names can be found below the corresponding photograph.All images were taken by individuals within the research team.

A PPE N D I X 3
A comprehensive list of all the species observed in the sample-based dataset, including the family to which the species belongs, their abundances at each of the three dive sites, and total abundance.The dive sites samples were Chains (C), Anton's (A), and Bikini (B).Unique species belonging to this dataset are indicated by an asterisk.

A PPE N D I X 4
A comprehensive list of all the species observed in the incidence-based dataset.The species names are displayed along with their respective families, the total number of individuals observed as well as at which of dive sites the observations were made.Unique species belonging to this dataset are indicated by an asterisk.

A PPE N D I X 5
A table summarising the three biological indices utilised in this study obtained for each sampling event at each dive site: The Shannon-Weaver index, the Simpson dominance index, and the uniformity index.The uniformity index is calculated for the average Shannon-Weaver value for each dive site.The dive site is denoted by the first letter of its name and the sampling event is denoted by a number next to the letter of the dive site.Nudibranch abundance as well as diversity is also depicted for each sampling event at each dive site.
over a period of 2 weeks during April and May 2022.Incidencebased (IB) observations from citizen scientists were collected across 8 dive sites on Two-Mile Reef and included the three sites for the SB observations (Figure 1, Appendix 1).The depth for the IB observations ranged from 11.80 to 32.60 m.Incidence-based observations from citizen scientists were collected from February to July 2022.

F
I G U R E 1 A map of Two-Mile Reef off Sodwana Bay, South Africa.The three study sites where nudibranch sample-based data were collected are Chains, Anton's, and Bikini (indicated in bold on the top right-hand corner).The dive sites where incidence-based observations were made are also listed on the top right-hand corner.Coloured circles on the map represent the dive sites listed.
of community 0.00 < E ≤ 0.50 Depressed 0.50 < E ≤ 0.75 Unstable 0.75 < E ≤ 1.00 Stable was determined to be 0.84, 0.86, and 0.87 for Chains, Anton's, and Bikini, respectively.The overall average Simpson dominance index for all three sites was 0.86.The average uniformity index for each dive site respectively was 0.64, 0.62, and 0.63 for Chains, Anton's, and Bikini.The average uniformity index across all three dive sites was 0.63.All data regarding the diversity indices are viewed in Appendix 5.
Venn diagram depicting the number of species shared by the two datasets as well as the number of species unique to the sample-based and incidence-based datasets.F I G U R E 4 Species accumulation curves with (a) generated from the sample-based dataset and (b) generated from the incidence-based dataset.variables that may influence sampling.This investigation could then be used to make recommendations for improved experimental design.Nudibranch communities studied within the SB dataset indicate high levels of diversity, as shown with the combined average Shannon-Weaver index value (Appendix 5).Dominance of species is portrayed by the Simpson dominance index which indicates that

F
I G U R E 5 Rank abundance curves with (a) generated from the sample-based dataset and (b) generated from the incidence-based dataset.TA B L E 4 A summary of the generalised linear model investigating the effect of environmental and experimental parameters (predictor variables) on nudibranch abundance (response variable).
such as those investigated by Goddard et al. in 2011 can result in ecological disruption and thus the importance of understanding how nudibranchs are influenced by environmental factors in order to predict future potential impacts on ecosystem functioning.
SB and IB datasets combined.Dominant species observed in the two datasets included H. wasinensis, G. bonwanga, and C. hamiltoni, which corresponds with other research indicating that these species are common along the east coast of Africa.However, there is a data gap for G. bonwanga as it is a newly described species and therefore life history data along with distribution for this species needs to be collected.The average pH value obtained in this study was in line with the global average and can thus be used as a baseline reading for comparison in future studies.Nudibranchs have the potential to be effective bioindicators of environmental changes as well as reef health.This has been shown through range shifts, along with regional declines in California and Australia.The environmental and experimental factors investigated here did not have an influence on nudibranch abundance, species richness, or composition.Therefore, future experimental procedures to quantify nudibranch populations should not prioritise the influence of the strength of the surge and current, depth, dive time, or number of observers.Future studies can improve their methodologies by increasing the temporal scale of sampling in order to accurately quantify communities and incorporate larger variation within environmental factors.These changes to the sampling procedures will enable long-term monitoring of nudibranch populations and allow for insight to be gained on the impact of climate change in Sodwana Bay, which may enable nudibranchs to be classed as an effective bioindicator organism.AUTHO R CO NTR I B UTI O N S L. Garner: Conceptualization (lead); data curation (lead); formal analysis (lead); investigation (lead); methodology (lead); software (lead); validation (equal); visualization (lead); writing -original draft (lead); writing -review and editing (equal).C. J. Oosthuizen: Conceptualization (supporting); data curation (supporting); formal analysis (supporting); funding acquisition (lead); investigation (supporting); methodology (supporting); project administration (lead); resources (lead); software (supporting); supervision (lead); validation (supporting); visualization (supporting); writing -original draft (supporting); writing -review and editing (supporting).